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Classical nucleation theory (CNT) is the most widely used framework to describe the early stage of first- 
order phase transitions. Unfortunately the different points of view adopted to derive it yield different kinetic 
equations for the probability density function, e.g. Zeldovich-Frenkel or Becker-Doring-Tunitskii equations. 
Starting from a phenomenological stochastic differential equation a unified equation is obtained in this work. 
In other words, CNT expressions are recovered by selecting one or another stochastic calculus. Moreover, it is 
shown that the unified CNT thus obtained produces the same Fokker-Planck equation as that from a recent 
update of CNT [J.F. Lutsko and M.A. Duran-Olivencia, J. Chem. Phys., 2013, 138 , 244908] when mass 
transport is governed by diffusion. Finally, we derive a general induction-time expression along with specific 
approximations of it to be used under different scenarios. In particular, when the mass-transport mechanism 
is governed by direct impingement, volume diffusion, surface diffusion or interface transfer. 


I. INTRODUCTION 


Over the last century the classical approach to nucle¬ 
ation has reported an excellent ability to catch the es¬ 
sential rules underlying noise-induced phase transitions. 
The set of ideas which constitutes this framework were 
developed for over more than half a century. This devel¬ 
opment started with Gibbs’ work 1 ® on near-equilibrium 
phase transitions. Several years later, |Volmer and We-] 
bei 1 ® introduced kinetic aspects improving the purely 
thermodynamical picture given by Gibbs. This was fur¬ 
ther developed from a more atomistic point of view by 
Farkas^ who developed Szilard’s ideas and by Becker and 
Doring^, Tunitskii®, Frenkel and Zeldovich 11 within 


the context of liquid-vapor transitions. Finally, Turnbull 
and Fisher extended such a formalism with the aim of 
describing solid nucleation firstly from liquid and then 
from solid phases. Besides being intuitively appealing, 
classical nucleation theory (CNT) has shown an over¬ 
whelming robustness. Not only is it able to generate a 
satisfactory estimation for the nuc leatio n rate equation 
which is good for practical purposed- 1 ® HI. but also it pro¬ 
vides a natural mechanism for cluster formation which 
turns out being more real than initially expectecf 1 ®!. In 
more recent years, a massive amount of studies have reno¬ 
vated, if not improved, CNT base d on e ither phenomeno¬ 
logical or fundamental grounds 7 111(7 33 . The primary 
goal of this paper is to describe a systematic method to 
derive the Fokker-Planck equation for the probability to 
find clusters of a given size within the general setting 
of CNT. This results into a surprising generalization of 
the usual form of the stationary distribution compared to 
what is usually assumed in the classical theory: rather 
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than a simple Boltzmann distribution depending on the 
work of formation of a cluster, there is a state-dependent 
prefactor that depends on the chosen stochastic calcu¬ 
lus. For the particular choice of the Stratonovich calcu¬ 
lus, the resulting expression reproduces the one recently 
derived from a more elabo rate formalism based on fluctu¬ 
ating hydrodynamic^^. Thus, one of the main contri¬ 
butions of this work is to provide a simpler route to that 
result based solely on CNT. From that, another impor¬ 
tant contribution comes up with the derivation of a gen¬ 
eral induction-time expression accompanied with specific 
approximations for it to be used in a variety of scenar¬ 
ios. Special attention is paid to those where the mass- 
transport mechanism is governed by direct impingement, 
volume diffusion, surface diffusion or interface transfer. 
To do this we followed the scheme outlined below. 


In section [TT] a stochastic differential equation (SDE) 
is proposed to model the time-evolution equation for the 
cluster size under a general stochastic calculus which is 
parametrized by a so that a. = 0, \ and 1 corresponds 
to Ito, Stratonovich and backward-ltd convention d 34 ! 35 !, 
respectively. The state-dependent components of the 
postulated SDE are derived from a phenomenological 
point of view, following the CNT reasoning. The Fokker- 
Planck equation connected to such a SDE is then intro¬ 
duced in section mu being now apparent the similari¬ 
ties with the Zeldovich-Frenkel equation, except for pre¬ 
senting an additional term which affects the effective en¬ 
ergy barrier. The effects of such a term are qualitatively 
explained in the same section. While the Zeldovich- 
Frenkel equation is recove red under backward-ltd con¬ 
vention, the one derived by LutskcJ ^1221 j s obtained under 
Stratonovich calculus. Lastly, the ltd convention is stud¬ 
ied which unveils the necessary condition for the different 


Fokker-Planck equations (FPEs) to converge. Section IV 


is devoted to the derivation of a generalized expression for 
the mean first-passage time as well as different approx¬ 
imations for the most common experimental conditions. 
Finally our results are summarized in section [V] 
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II. THEORY 


by means of the definition of the stochastic integral 


Nucleation is a quite complex process which certainly 
requires a huge amount of relevant variables (order pa¬ 
rameters) to be fully described. However, as we com¬ 
mented in section [TJ CNT is based on a manageable de¬ 
scription of the process which considers a single order pa¬ 
rameter, namely the size of the emerging embryo (cluster) 
of the new phase within the old one. This is commonly 
assumed to be spherical and to have the same physical 
properties as those of the final stable stat d 13 * 1 ^ . This ap¬ 
parently crude simplification is nevertheless more than 
efficient when it comes to predicting the functional de¬ 
pendence of the nucleation rate on the thermodynamic 
quantities involved 13 ^. That is why our study starts 
with a proposed SDE for the time evolution of the clus¬ 
ter size. The reasons why we are motivated to propose 
such an equation are mainly two. On the one hand, it 
is known that nucleation is a thermally activated pro¬ 
cess involving the escape fro m a metastable state via 
overcoming an energy barriei ^ 14 * 36 l Thus, if the nu¬ 
cleation process is governed by a single order parame¬ 
ter it should undergo, at least in good approximation, a 
Langevin equation as pro posed b y Kramers 37 . On the 
other hand, recent studied ^ 32 ! 33 ! have shown a formal 
derivation of such kind of SDEs when the mass transport 
is governed by diffusion. That said, a general Langevin 
equation is assumed and their drift and diffusion terms 
will be deduced from phenomenological arguments, and 
inspired by the ideas underlying CNT. Without further 
ado, we will proceed writing this up. 


A. Stochastic dynamics 

Let us consider a system which is set in a metastable 
state and direct our attention to an arbitrary spherical 
cluster of those randomly growing and shrinking. Let A' 
be the size of such a cluster, accounting for the number 
of molecules inside it. Following Kramers’ reasoning for 
thermally-activated escape processed 7 we will propose 
the model equation, 

dX(t) = r]- 1 (A, t)F( A, t)dt+ V 2k B Tr ] - 1 (X,t ) * dB(t) 

(1) 

with 77 being an effective viscosity, F being the effec¬ 
tive force acting on A, T is the temperature, ks the 
Boltzmann constant and dB(t) being a Wiener process. 
Intrinsic to this equation is the consideration of A' as 
a continuous variable, approach first introduced by |Zel- 
ldovicli l l " | l3 i 3i 3 ED. Based on the work of Ree et al. 13 ^ it was 
shown (e.g. Sec. 9.1 of Ref. EH that clusters of more 
than a few molecules can be fairly well described within 
the framework of the continuous approach, and so will we 
consider it as a good enough paradigm. The star product, 
*, was introduced to remark that we are using a general 
stochastic calculus (hereafter called a-calculus) defined 


I a [R(x,t)\ = f dB(t')R(x(t'),t') 

Jt 0 

n —1 

:=ms - lim V R{x(t*), t*)AB t ( 2 ) 

n—> 00 ' 

2=0 

where R(x,t) is a left-continuous function, i.e. a func¬ 
tion which is continuous from the left at all the points 
where it is defined, the symbol means definition, 

(ms-lim) represents mean squared limit, i.e. a second 
moment convergence, t* = (1 — a)ti + at i+ 1 and A B t = 
B(t i+ 1 ) — B(ti). It can be showiP^HH] that choosing 
a = 0, 4 or 1, one recovers Ito, Stratonovich or backward- 
Ito’s definition. At the moment we will focus our efforts 
on obtaining a good estimation of the drift term and 
later on we will discuss the consequences of selecting one 
or another value of a. 

As previously mentioned, our aim is providing an 
alternative route of derivation of CNT starting from 
the cluster-growth law. Based on CNT, we know that 
the effective force, F, is related to the work of cluster 
formation , W, through its derivative, 


F(X,t) 


dW(X,t) 

dX 




( 3 ) 


The work of cluster formation is usually expressed as the 
increment of free energy experienced by the system due 
to the emergence of a cluster of size A. Depending on 
the system under consideration that work is specified in 
terms of either the Gibbs (AC?) or Helmholtz (AA) free 
energy, or the Grand (Landau) potential (AH). In spite 
of making the derivation as general as possible, we do not 
specify a given thermodynamic potential since knowing 
the functional dependence of F on W is enough by far. 

The derivation of an expression for 77 requires nonethe¬ 
less a slightly longer discussion. According to CNT, the 
effective time a cluster will spend to lose a molecule is 
given by the inverse of the difference between monomer 
attachment, /, and detachment, g , frequencies, 


A- = 


1 

g(X,t) - f(X,ty 


( 4 ) 


Within nucleation regime, clusters experience a stronger 
force to shrink than to grow due to their metastable na¬ 
ture. This results in a higher detachment rate than the 
corresponding attachment frequency (Chap. 10 or Ref. 
[13} with tv -C Is, typically of order 10 ~ 7 — 10 _12 s as 
discussed by Kashchiev (Figs. 10.2, 10.4 and 10.6 of Ref. 

EH). 

On the other hand, from the definition of the Kramers- 
Moyal coefficient^^ 


D^(x,t) = Alim -([A (t + T)-x] n ) 
n\ t-s-o r 


( 5 ) 


X{t)=x 
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it was argued (Refs. HD and B7>|) that the first Kramers- 
Moyal coefficient related to equation (|T|) is given by, 

D^(X,t)=a(X,t) + a db ^b(X,t) (6) 

with a(X,t) = g~ 1 (X,t)F(X,t) and b(X,t) = 

a/ 2 A:bT? 7 _ 1 (X, t) being the drift and diffusion forces, and 
P = l/k B T. Now, we can use the fact that the time r<_ is 
expected to be small compared to the typical time scale 
associated with a significant change of (X(t)) so that 
the limit in Eq. § can be approximated by evaluation at 
r = 7 y_ for which (X(t + r<_) — X(t)} = — 1, giving 

D W (X, *) = ( ^ ) ~ t) - g(X, t ) (7) 


Using equation (10) into (12) we get, 

W'{X, t) 


f{X,t)-g{X,t)=-f{X,t)- 


k B T 


(13) 


d 


+ /(X, t) ^ In [cr(X, a; t)f(X, f)], 

where we replaced by “=” assuming that these ap¬ 
proximations are accurate enough for all practical pur¬ 
poses. Finally, by substituting equation (131 into 0 
and equating to equation 0 we eventually get, 


V 1 {X,t) = Pf(X,t), 
cr(X,a;t) = / Q_ 1 (X,f). 


(14) 


Having equations ([3]) and (14), the Langevin equation 
proposed at the very beginning (Eq. |T]) can be finally 
rewritten, 


which agrees with the cluster-growth law derived in 
CNTpH 


(^f) =f(X,t)-g(X,t) ( 8 ) 

\ az / CNT 

As we know from CNT, whereas the analytical e xpres- 
sion of /(X, t) can be derived from collision theorjff““l 
finding the frequency g(X 7 1) is not a trivial task. This is 
so because monomer detachment depends on parameters 
characterizing the cluster which may differ appreciably 
from those of the bulk new phase. To get rid off this 
problematic quantity we will follow the same reasoning 
as in CNT. 

Let us assume that f(X,t) is a well determined quan¬ 
tity. We can evaluate the difference f(X,t) — g(X,t ) in 
terms of /(X, t) making use of the detailed balance equa¬ 
tion, 


/(X - 1, t) P(X - 1, t) = g(X, t) P(X, t) (9) 

with P(X, t) being the quasi-equilibrium probability den¬ 
sity function^, 

P(X, t) = P 0 (t)a{X, a ; t) e - pw(X ’ t] (10) 

where Po(f) is an instantaneous normalization constant 
and er is some function of X, a and t. Looking to ac¬ 
complish our goal, we will follow Kashchiev 47 l and ap¬ 
proximate /(X— 1, t)P(X — 1, t) by the truncated Taylor 
expansion about point X, 

/(X-m)P(X-m) 

~/(X, t)p( X, t) - Jv/(X, t)P(X, t ), ( 11 ) 

and, hence, 

/(X, t) - g{X, t ) ~/(X, t) In /(X, f)P(X, t). (12) 


dX(t) = —/(X, t) dl3W g ^ : dt + y/2f(X, t) * dB(t). 

(15) 

This stochastic equation is equivalent to that derived by 
Lutskop 32 from fluctuating hydrodynamics for a single¬ 


order parameter, when we set a = | (i.e. Stratonovi ch’s 


calculus). Moreover, as it will be checked in section III 


this equation is statistically equivalent to the Zeldovich- 
Frenkel equation for a = 1 (i.e. backward-Ito’s calculus), 
given that it produces the same time-evolution equation 
for the probability density function. It is thus interesting 
to see how the hypotheses underlying Zeldovich’s deriva¬ 
tion are equivalent to choose a specific stochastic calcu¬ 
lus. Besides, now the equilibrium with the thermal bath 
is always ensured regardless the value of a, unlike the 
common belief that only the backward-Ito convention is 
capable of guaranteeing equilibriunP^. Actually, such a 
belief is reached after considering the equilibrium regime 
as equivalent to a Boltzmann distribution law. Never¬ 
theless from previous studies we know this is not true 
for nucleation, where the equilibrium distribution de¬ 
rived from a fluctuating-hydrodynamic framework shows 
a state-dependent prefactor ! 29 ! 3 ^ That is why we free the 
derivation herein presented of that restriction assuming 
instead a local-equilibrium law (Eq. 10 1 with a general 
state-dependent prefactor. This results in a general equi¬ 
librium distribution function (using Eqs. (10) and ( fl4| )) 
which yields a state-independent pre-exponential factor 
for a — 1, as expected. Noteworthy in such a case the 
resulting theory cannot be covariant as will be shown in 
section IIII Bl 

Nevert heless, equation (fl5| differs from that produced 
by Tunitskii’s equatiorP® 


dX(t) = ( f{X,t ) - g(X,t))dt+^/f{X,t) + g(X,t) dB(t) 

(16) 

which is interpreted under Ito’s convention, i.e. a = 0. In 
order to know under which circumstances both equations 
are equivalents, we will follow the same reasoning as that 
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led us to equation (14), and so we arrive at, 

f(X,t)-g(X,t)= -W'{X,t)r\X,t), 
f{X , t) + g(X, t) = 2k B Tr,~\X, t). 

Then, 


g(x,t) 

f(x,t) 


=i + 


l 




1 - 


( W'(x,t) \ 

\ 2fc b T ) 


k B T 


(17) 

(18) 


(19) 


To get the usual approximation of g{X , t) given in CNT 
(e.g. Eq. (10.90) of Ref. [T3]), i.e., 


we need 


9 (X,t) = f{X,t) exp 


d$W{X, t) 
dX 


(dpw{x,t) 

l, dX 


< 1 . 


( 20 ) 


( 21 ) 


given that in this limit one gets 

g(x,t) 


f(x,t) 


dpw(x,t) i ( dpw(x,t) 

X 2 \ dX 

( dj3W (X, t) 

^ dX 


' exp 


( 22 ) 


That way, our proposed model will recover the Tunit- 
skii equation under Ito’s convention for slowly-varying 
energy barriers, which indeed agrees with the hypotheses 
underlying CNT. 

Thus far, we have derived heuristically a model (Eq. 


15) which seems to be in accordance even with more 


rigorous and modern theories. In what follows, we will 
study the dynamics of the probability distribution func¬ 
tion (PDF) associated with equation (15) in order to 
make contact with CNT. This will make possible to get 
some important quantities such as the stationary distri¬ 
bution function (for undersaturated systems) or the nu- 
cleation (or induction) time. 


III. FOKKER-PLANCK EQUATION 


The time-evolution equation of the PDF, P(X,t), of 
the random variable X will be given by the the following 
Fokker-Planck equation®®! 


with 

Z(X,t) 


dP(X,t) dZ{X,t) 


dt 


dX 


(23) 


- <j f(X, t) — [pW(X, t) + (1 - a) In f(X, t)] 

+ f(X,t)^P(X,t). (24) 

- (f{X,t) m ^' t] +f(X,t)^j P(X,t), 

(25) 


with, 


mX, t) = /3W(X, t) + (1 - a) In f(X, t). (26) 

Now, the similarities between this FPE and that obtained 
in CNT are apparent. Yet more, the Zeldovich-Frenkel 
equation is recovered when the backward-Ito convention 
is adopted. Surprisingly enough, this naive model also 
reco vers t he FPE given in more recent rederivations of 
CNTtMH! when the Stratonovich calculus is considered. 


A. Short-time propagators: critical clusters with growing 
habits 


Now we are going to evaluate the impact of the extra 
logarithmic term. Since 0 < a < 1, it is evident that 
it entails an increase in the energy barrier with respect 
to the Zeldovich-Frenkel equation. However, there is an¬ 
other interesting effect coming from this additional term 
which has to do with the probability of a critical cluster 
(defined by j3W'(X*,t) = 0) to grow or shrink. While it 
is customary accepted that these probabilities must be 
the same, this is only true under Ito’s convention, as we 
will show right away. To this end, we will make use of 
short-time propagator^! 

It is known from collision theory that the analytical 
equation of /(X, t) is size-dependent and so is y/2f(X, t). 
Hence, their values will change from the initial to the fi¬ 
nal size during a unitary jump in the size axis, X. This 
change implies that the cluster feels different attachment 
rates in going from an initial size, Xo, to another, X. 
Here is where the choice of the stochastic calculus comes 
into play, since each of them corresponds to a different 
origin where evaluating the noise amplitude, which yields 
asymmetric probability distributions 43 44 for any a ^ 0. 
Considering a £ [0,1] the short-time propagato: 


re¬ 


lated to the FPE (22), and hence to the SDE (15), is 
given by the following equation @ 

p a (X,t + t|X 0 , t) 


( -ar dD(1 l{^ 


exp 


+a t 


ox 

2 d 2 P {2> (X a ,t) 


dX 2 


{X-X a -[D m {X a ,t)-2c,p( 2) ' {X a ,t)]r} 2 
4 rP( 2 )(X a ,t) 


2 JnTD(V(X a ,t) 


(27) 


with: 

D {1 \X a ,t) = — f(X a , t)(3W'(X a , t) + af'(X a ,t), 
DW(X a ,t) = f(X a ,t ), 

X Q = aX(f + r) + (1 — a)X(t) 

= aX + (1 - a)X 0 . (28) 


In the particular case of A' 0 = X* we know by definition 
that BW(X t ') = 0, i.e. the time evolution of the size for 
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ing the maximum deformation foi® a = 1, 

Paz=l{X,t + T\X*,t) 

{AX» + {f(,X,t)0W'(X,t)+f'(X,t)]T} 2 
_ 4T/(.Y,t) _ 

2v/t r/(A,t) 

(30) 

Given that /(A, t) is a monotonically increasing func¬ 
tion of X (e.g. Ch. 10 of Ref. IT5T> . it can be shown 
that such an asymmetry favours the growth, instead of 
the shrinking. This fact has been already reported by 
Brettschneider et al.® in a different scenario but with 
similar conditions. Inspired by their discussion, this re¬ 
sult has been schematically represented in Fig. [l] While 
this is a very interesting result which underlies the dif¬ 
ferent derivations of CNT, it is true that for very large 
critical clusters, A'* 1, the general short-time propa¬ 

gator becomes, 




p CNT (A ,t + T\X*,t) = lim p a (X,t + T\X*,t) 


exp -f (AX * )2 1 
P \ 4 f(X,,t)rJ 

2y/n f(X.,t) 


(31) 


FIG. 1. This figure is inspired by Fig. 4 of Brettschneider 
et al.® The top panel (a) is a schematic representation of 
the dependence of the monomer attachment frequency on the 
cluster size. Below, panels (b), (c) and (d) show qualitatively 
the effect of the logarithmic term on the short-time propa¬ 
gators, in particular when the initial condition is the critical 
cluster. It is readily observed how the symmetry of growing 
or shrinking is broken for all a ^ 0. 


a critical cluster is purely stochastic (Eq. [i~5| ). However, 
as we can observe in Fig. [T] the probability to grow only 
equals that to decrease for a = 0, given that 


where we have considered the Szilard limit 6 in order 
to approximate /(A, f) ~ /(A*,t). Thus, we see 
that the different a-calculi yield the same results only 
for very large critical clusters, which is equivalent to 
say for near-equilibrium systems. It is worth mention¬ 
ing that a recent rederivation of CNT from fluctuating 
hydrodynamicfpSEMU] (hereafter denoted dCNT) has for¬ 
mally shown that the right FPE is that given for a = |, 
otherwise this will not be a covariant theory. Hence, from 
the above discussion, the critical cluster will experience 
a slightly higher tendency to grow, unlike what would be 
initially expected from CNT. 


B. Equilibrium and stationary distributions 


p a =o(X,t + T\X*,t) 

- 1 PVT , ( ( A *») 2 ^ 

20r/(W*,£) eXP V 4 f(X*,t)Tj 


(29) 


where AA* = X — A*. This is in line with the short- 
time propagator directly derived from Tunitskii’s equa¬ 
tion (Eq. [l6| when /(A*,f) = g(X*,t), i.e. the probabil¬ 
ity to increase by one unit the cluster size must be equal 
to the probability to decrease by the same amount. 


In this section we are going to explore the solutions of 
the general FPE given by equation (221 and (25). It is 
widely known that finding the exact solution is a highly 
difficult problem, as well as potentially impossible. How¬ 
ever, as discussed by Hanggi et al.^1, an exact solution 
can be obtained assuming a stationary system with con¬ 
stant flux, 3 S , which is ensured by removing clusters once 
they reach a given size, A = X max > A* (see also Refs, 
ms an dEo}. Then, the steady-state distribution must sat¬ 
isfy P s (X max ) = 0. Let us consider f(X,t) = /(A') and 
<k(A, t) = C E > (A). Coming back to equation (25), 


Nevertheless, for all a £ (0,1] the short-time propaga¬ 
tor gradually becomes an asymmetric distribution reach- 


-vfs 



d/3<f>{ X) 
OX 


+ f ^8X 


P(A) 


(32) 
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which readily produces} 


r x 0 $(y) 

P S (X) =Ae~^ x ^ -a s e-^ A ')y ~J{Y) dY (33) 

with A being a normalization constant. When the bound¬ 
ary condition on P s (X max ) = 0 is imposed, Eq. (33) 
becomes, 


P,{X)=Xse~ p * (x) f 

■h 


X ma x „/3$(F) 


f(Y) 


-dY. 


(34) 


Now, we can note that for some change of variable X 
Y(X ), the stationary probability should fulfill, 


P S (Y) = P S (X )(g 


which imposes the following condition on f(X), 

!/(<*—i) 


/(F) = f(X(Y)) 


d_X\ 

dYJ 


(35) 


(36) 


Therefore, as we can check, the classical Zeldovich- 
Frenkel equation cannot be covariant 32 given the singu¬ 
larity occurring for a = 1. For all other values, whether 
the theory presents general covariance or not is condi¬ 
tioned by the definition of the attachment frequency. 
Thus far, the only derivation satisfying such a require¬ 
ment has been dCNT within the context of diffusion- 
control mass transport 32 . As results from such a study, 
the effects related to the non-covariant character of CNT 
are subdominant. 

If we consider an initial equilibrium (undersaturated) 
state, the stationary PDF becomes the equilibrium dis¬ 
tribution by identifying it with a zero flux regime. Then, 
imposing 3 s = 0 into equation (33) we arrive at, 


P eq {X) =Af( a - 1) {X)e- pw{ - x) = Ae~^ x] (37) 


The results introduced in this section will be remark¬ 
ably important when it comes to the derivation of the 
induction times, characterized by the mean first-passage 
time (MFPT). 


C. Semiadiabatic limit 

Nevertheless, real experiments could involve time- 
dependent coefficients. Following Weidlich and Haag^, 
we will assume that such a dependence of both / and >I> 
on time is controlled via a certain control function, n(t). 
Under this assumption we get, 

, df{X,K(t)) 

+ a OX ’ 

D^(X, K (t))=f(X, K (t)). (38) 


As evident, if k is time-independent we recover the pre¬ 
vious results. In most cases, the control function will be 
either the average density of the metastable state, tt(t) = 
Pav(t), or the temperature of the system, k(£) = T(t). 
In those cases in which «;(£) is a slowly-varying function 
of time, with respect to a typical time scale, one expects 
that the stationary solution follows adiabatically the mo¬ 
tion of k(£), i.e. the system reaches a stationary state 
almost instantaneously. This hypothesis is also known 
as semiadiabatic limilP 1 . The time scale that character¬ 
izes this limit can be interpreted as the relaxation time 
toward the initial metastable state, as pointed out by 
Talkner and Luczka^ 52 l Under these circumstances the 
zeroth order approximation for the quasi-stationary PDF 
is 


Pqe(X\t) ~A(£) e -P*(Xrt*)) 
_ 3 s e-m x Mt)) 



e P<S>(Y, K (t)) 

f(Y,K(t)) dY ’ 


(39) 


and so, the quasi-equilibrium PDF for undersaturated 
conditions, for which the flux necessarily vanishes (3s = 
0), is 


P qe {X-,t) ~A(t)e-P* {x ’ K W\ (40) 


The semiadiabatic approach is quite useful in order to 
make a first approximation of the MFPT under non¬ 
stationary conditions. Nonetheless, for a more accurate 
approximation of induction times the path integral for¬ 
malism developed by Getfert and Reimann 53 -l should be 
considered. 


IV. ESTIMATION OF INDUCTION TIMES AND 
NUCLEATION RATES 

In this section an approximation for the MFPT and 
so for the nucleation rate will be provided for most typ¬ 
ical mechanisms governing mass transport, i.e. the at¬ 
tachment rate. To do that, we will firstly get an ap¬ 
proximation of the MFPT and then we will particularize 
it by using the most used expressions for f{X) in the 
literature 13 . The reason why we will focus on obtaining 
an approximation for the MFPT (hereafter denoted by 
t) is that this quantity is closely linked to the nucleation 
rate, J. Indeed, within the range of applicability of CNT 
one expects (e.g., Sec. 10.6 and Eq. 72 of Ref. l54l) . 

da. ^ Pav^~a * (dl) 

where the subscript a has been introduced to highlight 
the fact that it will depend on the a-calculus selected. 
That way, r can be understood as the mean time required 
for nucleation to occur. For this purpose we will follow 
Risken^ to get an approximation of r. To begin with, 
we will consider stationary conditions, i.e. / and W time- 
independent. From the results thus obtained, a naive ap¬ 
proximation for time-dependent conditions will be read¬ 
ily produced by considering the semiadiabatic limit. It is 
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easy to show that equation (32) can be rewritten as, 


-Js — 


/( X)e~ mx) -^ (e mx) P(X) 
a { 


= _ e -P<p(x) 


0 /3a> 


dX 


with 


<p(X) = W(X) - ak B T\n(f(X)). 


(42) 


(43) 


If the barrier is relatively high, £5 S is expected to be very 
small. Hence, we can integrate equation (42) from the 
minimum, X = X min , to the maximum, X = X max , size 
allowed for a cluster, 


r X„ 


lx... 


e^ s) ds 

= P(X min ,t)e^ x ” 
~ P{X min )e^ x — 


1 - 


P(X. 


max ) 


P(Xmin) 


(44) 


Under these conditions, as (f s is assumed to be very 
small we can use Eq. (33) to approximate P(X m i n ) ~ 
and p(X) ~ That way, the dis¬ 

tribution function near X rnln will be approximately given 
by (see Ref. 


(45) 


P(X) ~ P{X mm )e-^ x) -^ x ^l 


Then, we can get the following expression for the 

MFTFfSH 


-(a) 


P(X min ,t)e^ x ^ /*™“ e~ ms) ds 
P(X min ,t)eP*( x ^)/ f* ma * ePvWds 


nXrr 


>Xrr 


-ms) ds 


r x „ 


/ x„ 


e Ms) ds. (46) 


The usual procedure 34 ! is to expand both exponents 
around X m i n = 0 and X* respectively, but in this prob¬ 
lem such a method cannot be used, since /(X) goes to 
zero as X —► d 13 l 14 ( Nonetheless, the main contribution 
to the first integral stems from the region around X rn ,j jTl . 
Thus, in the case of 3D nucleation the free energy term 
will be governed by the surface term 13 and hence, the 
exponent involving /34> can be truncated in good approx¬ 
imation as below 

P$ 3 d(X) = p (ex 2 / 3 + (1 - a)fc s Tlog(/(X))) , (47) 


where 0 would be a measure of the surface tensiorP^ HI 
On the other hand, for 2D nucleation, the work of clus¬ 
ter formation near X min is governed by the line-tension 
term, 

/3$2d(X) = p (OX 1 / 2 + (1 - a)fc s Tlog(/(X))) , (48) 


with 0 being here the scaled line tension. Hence, we get 

fX-max pOO _ 

/ e~W a) ds~ / e- p *Wds, (49) 
Jx min Jo 

whose exact results are collected in table |T] Note that 
these results involve the Euler gamma function, T(n). 
Besides, 

p Xmax pOO 

/ e^ s) ds ~ (s-x*) 2 } ds 

J x,„,,„, Jo 


1 


:X e ‘ 




with 


i _ I pW'{x.)\ 
A V 2t r 

X =A f 1 + erf 


u^ + W"(X*)\ 

2tt ksT 




(50) 

(51) 

(52) 


The pre-factor A~ 4 turns into the well-known Zeldovich 
factor,^ Zd, when one selects a = 0, i.e. Ito calculus. In 
fact x _1 could be considered as a generalization of the 
Zeldovich factor since it plays the same role as the latter 
in the classical expression of MFPT^ 


7(*. 


(53) 


The integral of equation (49) depends on the monomer 


attachment mechanism. In order to perform such an in¬ 
tegral we have used the expressions given in the literature 
for the most usual experimental situations (see table [I]) 
for both Homogeneous (HON) and Heterogeneous nucle¬ 
ation (HEN): 

• f(X) = CX 1 / 2 for 2D HEN of clusters with mono- 
layer height, with £ given by equations (10.6), 
(10.7), (10.63) and (10.66) of Kashchiev^ 

• f(X) = £X 4 / 3 for 3D HON and HEN of caps in 
liquid or solid solutions, with £ given by equations 
(10.18), (10.20) and (10.24) of KashchievS® 

• /(X) = £X 2 / 3 for 3D HON and HEN of caps when 
the monomer attachment frequency is controlled by 
direct impingement or by interface transfer, with 
£ given by equations (10.3), (10.4), (10.5), (10.9), 
(10.60), (10.61), (10.64) and (10.65) of Kashchiev 13 


Hereafter the integral of equation ( [49] ) will be called 1(a). 
From equation (46) we can finally give the approximation 
of the MFPT, 

r(«)~l X 7(a)/(X s )-V^, (54) 

which makes it possible the comparison with that pre¬ 
dicted by CNT (Eq. [54]), 


r(a) _ z^jl + aO ((X t )-*)]-l f(X t )-<*I(a) 

T CNT Z d f(X*) 1 

=I(a) /(X*) 1- “[l +aO((X 4 )- 2 )]-), (55) 






















which under backward-Ito calculus turns into 


r(a = 1 ) 


1 ( 1 ). 


(56) 


The fact that the approximations are not exactly the 
same, even though a = 1 produces the same FPE, it 
is because we followed a different route to derive r. 

One immediately observes the similitude of equation 


(54) with Kramers’ law. In fact, the former is in accor¬ 


dance with the latter with a pre-factor deduced analyti¬ 
cally. Yet more, equation (54) allows to reach an approx¬ 


imation of the nucleation rate via equation (41). Now 


this result can be extended as a first-order approxima¬ 
tion for time-varying conditions, as shown by Getfert and 
Reimann^ for slowly time-dependent Kramers-Moyal 
coefficients. Thus, under the assumption of the semia- 
diabatic limit, and following the notation introduced by 
Getfert and Reimann^, the instantaneous MFPT can 
be estimated as follows, 


r(a;t) ~ ± x (t) I(a, «(*)) f(X., K (t))~^ x ^\ 


(57) 


While these results can be very useful in order to char¬ 
acterize the time required to start the phase transition, 
we must bear in mind they are approximations. To get 
more accurate estimates of this characteristic time an ex¬ 
act numerical integration of the above equations can be 
performed. Indeed, the best estimation will be deter¬ 


mined via stochastic integration of equation (151. The 


main advantage of a numerical approach is that it is 
valid both for time-independent and time-dependent co¬ 
efficients. 


within the range of applicability of CNT, since the initial 
state is near equilibrium. Besides that achievement, we 
employed the tools of the theory of stochastic processes 
to go from the stochastic equation to a unified FPE. It 
is called “unified” given the fact that it contains both 
the classical Zeldovich-Frenkel equation when we select 
backward-Ito convention, and the one derived by Lut- 
sko and Duran-Olivencia 32), for the usual Stratonovich 
calculus. 

Another interesting result generated by this stochas¬ 
tic formalism was found while studying the short-time 
propagators. They constitute the tools required to know 
the probability for a cluster to grow or shrink. Surpris¬ 
ingly, when we set the initial cluster to be the critical 
(whose growth law is supposed to be purely random), 
one finds that all the interpretations of the noise will 
give a tendency to grow. It could be argued that this 
fact would be enough to make us select Ito interpreta¬ 
tion. Nevertheless, the only interpretation which has re¬ 
ported the ability to produce a covariant theory has been 
the Stratonovich calculus^. Yet more, the study on the 
general covariance of the theory reported another inter¬ 
esting result, namely Zeldovich-Frenkel equation cannot 
be fixed to be covariant. This result indeed sheds some 
light on the question whether or not it is worth trying to 
fix CNT by considering much more sophisticated models 
for cluster. 

Finally, estimates of the nucleation time and rate were 
computed. The approximation we reached for these 
quantities were specifically applied to the most usual ex¬ 
perimental situations. An inevitable similarity to the 
CNT expressions appeared. Lastly, a comparison with 
such an expression was performed by computing their 
ratios. 


V. CONCLUSIONS 
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Appendix A: Short-time propagator 

In order to get the expression of the short-time propa¬ 
gators in a general a convention, we will follow the proce¬ 
dure introduced by Wissel^, which can also be consulted 
in RiskerM 

The formal solution of a general FPE can be written 
as (e.g. Ref. flTTli 


p(x, t\x ', t') = T exp 



£fp[*(s), s]ds 


8(x 


x'), (Al) 
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TABLE I. Integration of equation (491 by employing some of the most used expressions for the monomer attachment rate and 
the corresponding approximations of the MFPT. 


Nucleation type mass-transport mechanism f(X) 


f™e~^ s) ds 


LvIFPT 


2D HEN 

disks with monolayer height 


direct-impingemenlT 

interface-transfeJ b l^ 


Cl* 


1/2 


h(a) 


_ 2r( a +i)Ci 

— RO + I QO + 1 




3D HON/HEN 

spheres/caps from 
liquid or solid 


volume-diffusio: 


iofl C 2 * 1/3 Ha) = 

■— 1 2/3 2 © 2 


f/2 X e^ (Jf * ) 


3D HON/HEN 

spheres/caps from vapor, 
liquid or solid 


direct-impingementr 
interface-transfeirT^ 


C 3 a 2/3 


r ( \ - 3r(?^±i) c 3 “ _1 

/ 3 (a) — 2 a + 1 2 a + 1 


2(3 2 © 2 


\h X e Pv(x ^ 


3D HON/HEN 


arbitrary 


(X s/v 


I(a) - 3 r ( l (^-§+ 1 ))^~ 1 l Iye ^( x .) 


a Cl prefactor multiplying N 1 / 2 in Eqs. (10.6) and (10.7) of Ref. [13] 

b Cl prefactor multiplying N 1 / 2 in Eqs. (10.63) and (10.66) of Ref. [HI 

c C 2 prefactor multiplying /V 1 / 3 in Eqs. (10.18), (10.20) and (10.24) of Ref. [HO 

d C3 prefactor multiplying N 2 / 3 in Eqs. (10.3)-(10.5) and (10.9) of Ref. [13] 

e C3 prefactor multiplying N 2 / 3 in Eqs. (10.60)-(10.65) of Ref. [13.: 


where £pp denotes the Fokker-Planck operator (Eq. A3) 
and T is the time-ordering operator. For small time in¬ 
tervals, t = t — t', equation (Al) can be approximated 
by 

p(x,t + t\x', t ) 

= [l + C FF [x' + aAx, t\r + C9(r 2 )] S(x — x'), (A2) 

where we have used the stochastic integral"^ introduced 
in equation ([2|. Now we can perform the differentiation 
in the Fokker-Planck operator 

a o2 

£fp (u(x),t) = D [ 1 ) [u(x)+\ + ^-^Z? ( 2 ) [u(a;),t], 

(A3) 

with u{x ) = ax + (1 — a)x'. Accordingly one obtains 

gn „g 2 r)( 2 ) 

C FP (u(x),t) = - a ^ ( u(x),t ) + a 2 ( u(x),t ) 


dx 2 


pi n(2) 

D^tidx), t) — 2a —-—( u(x),t ) 
ox 


+ D (2) {u{x),t)-^. 


d_ 

dx 

(A4) 


By substituting equation (A4) into (A2), and after that 
using the Taylor expansion of the exponential function 
we get, 

p{x,t + T\x',t) =e CFFMx) ’ t)T 6(x - x'). (A5) 


With the aim of obtaining an exponential function, the 
Fourier representation of the Dirac-i5 function will be 
used. Thus, we finally obtain the sort-time propagator: 


p(x,t + t\x' ,t) = 

( -ar afl ‘ l y i) ’ t) \ 

exp 

{x-x'-[D < ' 1 \u{x),t)-2a^XL(u{x),t)] T } 2 
\ 4 tD( 2 ) (u(x),t) / 


y / 47rrD(2T(^^~/y 


(A6) 


As can be checked, this short-time propagator recovers 
those presented by Dekker 1 '^ and Risken^ (e.g. Eqs. 
(4.55) and (4.55a) of Ref. [34]) for a = 0 and a = 1, 
respectively. This derivation can be seen as an alternative 
route to that carried out by Lau and Lubensky^l. 
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